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We study the elastic properties of phantom networks of Gaussian and nearly-Gaussian springs. 
We show that the stress tensor of a Gaussian network coincides with the conductivity tensor of 
an equivalent resistor network, while its elastic constants vanish. We use a perturbation theory to 
analyze the elastic behavior of networks of slightly non-Gaussian springs. We show that the elastic 
constants of phantom percolation networks of nearly-Gaussian springs have a power low dependence 
on the distance of the system from the percolation threshold, and derive bounds on the exponents. 
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I. INTRODUCTION 



Rubber and gels are large polymeric solid networks 
formed when polymers or monomers in fluid solutions are 
randomly cross- linked by permanent bonds. This process 
is called vulcanization or gelation, when the latter term 
usually applies to cross-linking of monomers or very short 
polymers — gels; while the former term usually describes 
formation of dense networks of long polymers — rubber. 
Rubber and gels are much more flexible than ordinary 
crystalline solids and, moreover, may remain in the linear 
elastic regime even in response to deformations increasing 
their dimensions far beyond their original, unstrained, 
size. Such a behavior is attributed to the network struc- 
ture of these materials, and to the fact that the elastic 
restoring forces are of entropic, rather than energetic, ori- 
gin. The simplest theory of rubber elasticity which cap- 
tures these essential physical features, is the "phantom 
network" (PN) model introduced by James and Guth . 
This model assumes that the configurations of the differ- 
ent polymer chains are independent of each other, and 
neglects the excluded volume interactions between the 
monomers. With these simplifying assumptions one can 
treat each polymer chain in the network as an ideal one. 
By averaging over the positions of the monomers one 
finds that the probability density of finding chain ends 
separated by f, takes a Gaussian form ~ exp [— \Br 21 \, 
where B usually depends on the temperature T. The 
free energy of the chain is proportional to (minus) the 
logarithm of this probability density and, therefore, pro- 
portional to r 2 , as if it is a linear spring of vanishing 
equilibrium length, which will be called Gaussian spring. 
In the PN model, the thermal averages of some quanti- 
ties can be calculated analytically due to the Gaussian 
form of the statistical weights and this makes it an 
excellent starting point for models with excluded volume 
interactions and entanglements 

The problem of gel elasticity introduces an additional 
complication already at the level of the PN model. In 
gels the network strands are very short and do not neces- 
sarily resemble Gaussian springs. Nevertheless, one may 
still construct a Gaussian model of gel elasticity, simply 
by replacing each bond of the gel by a Gaussian spring. 



In the absence of excluded volume interactions, the va- 
lidity of this model is justified by the fact that even if 
the elementary pair potential between bonded atoms is 
very different from that of a Gaussian spring, the effec- 
tive interaction between somewhat more distant atoms is, 
almost always, quadratic. This is a well known feature 
of long polymer chains [|| , but it has also been demon- 
strated for more complicated networks ||. De Gennes 
used an analogy between elasticity of networks of Gaus- 
sian springs and conductivity of random resistor networks 
and argued that rigidity, just like conductivity, ap- 
pears at the connectivity threshold, when a macroscopi- 
cally large network spans the system. He further argued 
that at the phase transition the shear modulus and the 
conductivity should have the same dependence on the 
distance of the system from the connectivity threshold. 
Surprisingly, the details of the argument of dc Gennes 
have never been worked out, i.e., there is no detailed 
calculation of the quantities characterizing the elastic 
response of Gaussian networks, namely the stress and 
elastic constants tensors. [There are several analytical 
studies of the statistical properties (including the elastic 
properties) of systems of Gaussian springs , but none 
of them makes such an explicit calculation.] In section 
H of this paper we derive exact results for the stress and 
elastic constants of Gaussian networks. We prove that 
the stress tensor of a Gaussian elastic network is equal 
to the conductivity tensor of an equivalent resistor net- 
work. A detailed proof of this equality, which holds for 
a Gaussian network of arbitrary topology, is given in the 
appendix of the paper. We also show that the elastic con- 
stants of a system consisting of a single spanning cluster 
of Gaussian springs vanish. We discuss the effect of the 
finite clusters which model the small molecules formed 
in the process of crosslinking and show that they play a 
crucial role in stabilizing the system. 

investigate the elastic behavior of 



In section III 
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phantom networks of nearly-Gaussian springs, whose en- 
ergy dependence on their extension includes a small quar- 
tic term additional to the quadratic one. A perturbative 
analysis yields an expression for the elastic constants. In 
section [V we use this expression to evaluate the elastic 
constants of phantom percolation networks ||, close to 
the percolation threshold p c . We conjecture a universal 
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scaling law for the elastic constants, C ~ (p — p c ) 9 , and 
derive exact bounds for the scaling exponent g. Section 
includes a short summary and discussion of the main 
results. 



II. ELASTICITY OF SYSTEMS OF GAUSSIAN 
SPRINGS— EXACT RESULTS 

A. Definitions in the theory of elasticity 

The theory of elasticity describes deformations of ther- 
modynamic systems in response to external forces. At 
a finite temperature, it is convenient to consider homo- 
geneous deformations of the boundaries of the system, 
which can be described by a constant matrix, My . When 
the system is strained, the separation between a pair of 
surface points, which prior to the deformation was R, 
changes to 

n = MijRj, (i) 

where the subscripts denote Cartesian coordinates, and 
summation over repeated indices is implied. Usually the 
energy of the system depends on the relative distances be- 
tween the atoms. The squared distance in the deformed 
system is equal to 

r 2 = r k r k = \l,,,n,M kl Rj 

= (M*M) y RiRj = (S l3 + 2rfe) R i R j , (2) 

where M* is the transpose of M, and rjij is the strain ten- 
sor, while Sij is the Kronecker delta. The strain tensor 
vanishes at the undeformed reference state. Expanding 
the mean free energy density in the strain variables 

f({v}) = /({0}) + PijVij + -^CijklVij'nM + ■■■ , (3) 

we identify the coefficients cry as the components of the 
stress tensor, while Cyy are the elastic constants (some- 
times referred to as the elastic stiffness tensor). 

The elastic constants of a thermodynamic system are 
related to each other through certain equalities. The ac- 
tual number of independent elastic constants depends on 
the symmetries of the system. Isotropic systems, for in- 
stance, have only three different non-vanishing elastic 

Constants: Cxxxx — ^yyyy ~ C zzzz — Cn ; Cxxyy — 
Cyyzz — Czzxx — ■ ■ ■ — Cl2 ) and C X yxy Cyzyz — 

C Z xzx = ... = C44 . Moreover, these three elastic con- 
stants obey an additional relation C\\ — C12 = 2C44, 
which reduces the number of independent elastic con- 
stants of isotropic systems to two. Frequently, one finds 
it more useful to describe the elastic behavior in such 
systems in terms of the shear modulus, /1, and the bulk 
modulus, k, defined by jl0| 

fi = C 44 - P, (4) 



and 

— / l^n + ^12), f° r 2 — dimensional systems 

\ f (C11 + 2C12 + P), for 3 — dimensional systems 

(5) 

where P = —tJ xx = ~ °~yy = ~ a zz is the pressure. When 
k and (i are positive, the system is mechanically stable 

B. Description of the system 

We consider a d-dimensional system shown schemat- 
ically in Fig. [p. The black circles in Fig. |] represent 
atoms while the zigzag lines indicate the bonds, attrac- 
tive pair-potentials, which connect them in a certain fixed 
(quenched) topology. Atoms which are found inside the 
volume of the systems are called internal atoms. Surface 
atoms have fixed coordinates on the boundaries of the 
system. The bonds connect atoms into clusters. Clus- 
ters containing only internal atoms are free to move in 
the entire volume. Cluster with both internal and sur- 
face atoms are non-free. Among them, one (and, in some 
cases, several) may extend from one side of the system 
to the opposite side. This is the "spanning" cluster. 




FIG. 1. A schematic picture of a network of springs. The 
system includes a spanning elastic network as well as some 
finite clusters. Atoms can be either internal, i.e., free to move 
inside the volume, or external, i.e., attached to a permanent 
positions on the boundaries. Non-free clusters have at least 
one external atom. 

The system which we study in this section consists of 
point-like atoms connected by Gaussian springs. The en- 
ergy of each Gaussian spring is given by 

(R a - i?' 3 ) = l -K a(i (R a -R^y ' = X -K a(i (R^f , 

(6) 
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where R a and R@ denote the positions of atoms a and /?, 
and R a ® is the distance between these atoms. The spring 
constant K af3 is assumed to have a fixed, temperature- 
independent, value. The total elastic energy is given by 
the sum over the energies of all the springs 



E 



(a/3) 



i/3> 



C. Elasticity of the system 

The components of the stress tensor of our system are 
related to the pair-potentials, <f> a /s(R a ^), via relation 




NkTSi 



(7) 



which has been derived thirty years ago by Squire, Holt 
and Hoover jL2| as an extension of the Born and Huang 
theory of elasticity jl3| , to systems at finite-temperature. 
In expression (f^) summation over all distinct pairs of 
atoms, a/3, is performed, where R"^ and i?"' 3 are the ith 

and the jth Cartesian components of i? Q/3 = R a — RP . 
The symbol ( } indicates a thermal average, while N and 
V denote the number of internal atoms and the volume of 
the system, respectively. For potential (^) the expression 
(^) reduces to 



0~iA = — 



K a ?RfRf 



v(a/3> 



NkT 



(8) 



where the sum is over the connected pairs. 

The two terms in the expression (||) are called the con- 
figurational and kinetic terms, respectively. The config- 
urational term can be divided into terms, each one in- 
cluding the sum over the bonds of one distinct cluster. 
Since there are no excluded volume interactions, these 
terms are independent of each other (the clusters do not 
interact with each other), and the contributions of the 
different clusters to the stress are additive. We identify 
the stress applied by each cluster as 
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N[kT 



\ (a (3) E cluster 



v 5ij ' 



(9) 



where 2Vj is the number of internal atoms of the cluster. 



free cluster (fc) should be as of a point-like atom. To 
prove this result (which is general and does not depend 
on the particular form of the pair-potential), we use the 
fact that for a free cluster (fc), one can integrate out 
d degrees of freedom (of say, R l ) in Eq.(||), and ex- 
press the terms appearing in it in the relative coordinates 
Rf = Rf - R\ {a = 2,..., Ni}. (This statement is cor- 
rect only in the thermodynamic limit, when the linear 
size of the system becomes much larger than the radius 
of gyration of the free cluster.) One can easily verify that 
in the relative coordinates Eq.(|^) may also be written in 
the following way 



\a=2 



8E 



NjkT 
V~ 



which from the equipartition theorem gives cry = 
— (kT/V)Sij. The stress applied by all the free clusters 
is simply 



N kT t 

v 1 



(10) 



where Nq is the total number of free clusters. Similarly, 
the contribution of the free clusters to the elastic con- 
stants is also as of an ideal gas, given by the kinetic term 

C3 



fc _ 2N kT 

°ijfcl - y d U°jk 



(11) 



E. Elasticity of the spanning cluster 

The stress and elastic constants of the spanning net- 
work of Gaussian springs with temperature-independent 
force constants, are temperature-independent. The free 
energy F of the spanning network is a function of the 
temperature T and the positions of the surface atoms 
{R s }. If the values of these variables change quasi- 
statically, then 



dF = -SdT 



£ ^£ct ' dR s 



(12) 



where S is the entropy, /,f xt is the external force which 
drags the surface atom s and summation is made over 
all the surface atoms. In a quasi-static process, the force 
/e Xt is balanced by the force f s applied by the network 
on atom s, namely 



D. The contribution of the free clusters 

The gas of free clusters is an ideal gas. Since the 
different clusters do not "feel" each other, it is intu- 
itively clear that the contribution to the stress of each 



fLt = f s = (J2 Kas 




(13) 



where summation is over all atoms a connected to atom 
s. The terms appearing in the thermal average in Eq.dlJ 
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are linear in the coordinates R a . Since the Boltzmann 
weight is a Gaussian, i.e., an exponent of a quadratic 
form of the coordinates, these averages coincide with the 
most probable values, namely their values at the ener- 
getic ground state, and therefore do not depend on the 
temperature. We thus conclude that f s is a temperature- 
independent quantity, and from Eqs.jl^) and ( |l3| ) we 
readily find that 



8 2 F 
8TdR s 



' dT 



0. 



The last result implies that F can be decomposed into 
two parts 

F(T,{R S }) = F 1 (T) + F 2 ({R S }). 

If we consider homogeneous deformations we may define 
a reference system and use the strain variables in- 
stead of {R s } 

F = F 1 (T)+F 2 ({r Hj }), 

The stress and elastic constants are the coefficients in the 
{?7}-expansion of F 2 [see Eq.(|])]. Therefore, they do not 
depend on the temperature. 

The stress applied by the spanning network is equal to 
the conductivity of a resistor network with the same topol- 
ogy. The stress of the spanning cluster (spc) [Eq.(ft)] 



spc 




NikT 

V tj ' 



can be rewritten in the form 



a—1 J , 




2^ K as RlR° a 

\ (as) £ spc 



(14) 



where the first sum is over all the internal atoms, while 
the second sum is over all the bonds connecting inter- 
nal and surface atoms. (The subscripts s and a denote 
surface and internal atoms, respectively.) In the thermo- 
dynamic limit we deduce from the equipartition theorem 
that the first and the third (kinetic) terms in Eq.(|l4|) 
cancel each other. We are thus left only with the second 
term 
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V 



K as Ri ) 



(as) E spc 



(15) 



The thermal averages in Eq. (|l5|) are of quantities which 
are linear in the coordinates of the internal atoms and 
therefore may be replaced by the equilibrium values of 
these quantities (see earlier in this section) . The equilib- 
rium values of R a minimize the energy of the spanning 
cluster 



E* 



\x a(i (R a - R f 



(a/3) £ spc 
d 



E 



E l^(Rf-R^ 



(a/3) G spc 



(16) 



The dependence of E spc on the components R" corre- 
sponding to one Cartesian direction, j, is included in the 
term E spc ■ The problem of finding the equilibrium val- 
ues of R a decouples into d scalar problems of finding the 
equilibrium values of R" . In order to calculate these val- 
ues we need to solve d sets of linear equation (one set for 
each Cartesian component): 







K 



a/3 



R«-RP 



= 0, 



(17) 



corresponding to the vanishing of the jth component of 
the force acting on each internal atom. (For each atom 
a, summation in the relevant equation is over all atoms 
(3 connected to it.) 

Let us define a resistor network with the same con- 
nectivity as the elastic network, in which each spring is 
replaced by a resistor with conductance K a<3 . The values 
of the electric potential at the internal nodes, {<p Q }, are 
obtained by minimization of the heat power produced in 
the network, P — J2( a /3) K af} {ip a - ip f3 ) 2 . Except for a 
prefactor of |, P is identical with Ej Spc (|l6|), where ip a 
plays the role of i?". If we replace R" by ip a in the force 
equations ([l?]), we obtain the set of Kirchoff equations 
enforcing the vanishing of the sum of currents entering 
the internal nodes of the network. By replacing R" by 
ip a , we define a mapping of the mechanical problem to an 
electrostatic one. In fact, we have d different electrostatic 
problems corresponding to each Cartesian component of 
the mechanical problem. They differ from each other 
in their boundary conditions, namely the values of the 
electric potential on the surface nodes, {<p s }. In the jth 
electrostatic problem, we set ip s equal to i?*, i.e., we as- 
sume that the electric potential at each boundary point 
is equal to the jth Cartesian coordinate of the point. 

The interesting question now is what is the analog of 
the stress tensor in the electrostatic problem. This ap- 
pears to be the conductivity tensor, Ey defined by 

(ji) = T.ij(Ej), 

where (j) and (E) are the volume averages of the cur- 
rent density and the electric field, respectively. More 
precisely, if we follow the mapping defined above we have 
the exact equality 
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(18) 



F. The stability of systems of Gaussian springs 



A detailed proof of this equality is given in the appendix 
to this paper. Here we just note that the proof consists 
of two steps: In the first step we show that in the jth 
electrostatic problem, because of the choice of boundary 
conditions, (E) is a unity electric field pointing in the 
(— j)th direction. In the presence of such an electric field 



(ji) 



-Ejj. On the next step of the proof we show that 



(ji), and therefore Ey, are given by the electrostatic 



equivalent of Eq. fll5| ) 



E 



(as) £ spc 



y^SpC 1 - 

ij _ y 



and therefore Eq.(|l8|) is valid. 

The elastic constants of the 
ish. We have already shown that C. 



(19) 



spanning network van- 



ijkl j 



the elas- 



tic constants of the spanning cluster of Gaussian 
springs with temperature- independent force constants, 
are temperature-independent. Therefore, we may calcu- 
late them at any temperature, and in particular at T = 0. 
At zero temperature the free energy coincides with the 



internal energy, given by Eq.(16|), where {R a }, the posi- 
tions of the internal nodes, take their equilibrium values. 
Suppose now that the system is homogeneously strained. 
The positions of the surface nodes, {R s }, change accord- 
ing to the linear transformation ([!]), with a constant ma- 
trix Mij. [Transformation (|l|) was originally defined for 
the separation between surface points. However, we can 
always set the origin of axes to be on the original (un- 
strained) surface, and in this case the transformation 
applies to the positions of the surface points.] In or- 
der to find the new equilibrium positions of the internal 
atoms, in the strained system, we need to solve the set of 
equation (|l7| ) with the new boundary conditions. Since 
both the equations and the transformation of the bound- 
ary conditions are linear, the new solution is given by 
rf = MijR". The elastic energy of the strained span- 
ning cluster is given by [see Eqs.(|) and ©] 



- J2 K^tf*? 



(as) £ spc 



K afj |~(M* M)ij Rf^RjP 



{as) £ spc 



£ K^\(2 Vl3 +8^RfRf 



(cxs) £ spc 

This gives the dependence of E on the strain variables, 
which include only linear terms in rjij. Since the elas- 
tic constants are the coefficients of the quadratic terms 
in the {^-expansion of the free energy [Eq.([j|)], we con- 
clude that 



(20) 



We have mentioned earlier in this section that sta- 
ble solid thermodynamic systems have positive bulk and 
shear moduli, k and fi [Eqs.(^) and (||)]. In phantom 
systems, the contributions of the spanning cluster and 
the ensemble of free clusters to k and fi are additive. 
Due to the vanishing of the elastic constants of the span- 
ning cluster (EQ), we find that its contribution to the 



,ps P c > o and 



(two- 



elastic moduli is: /i spc 
dimensions) or K spc = P spc /3 < (three-dimensions) 
[P spc is the negative (stretching) pressure applied by the 
spanning cluster]. The fact that k is not positive means 
that the spanning cluster alone is not stable against ho- 
mogeneous volume fluctuations. The contribution of the 
free clusters to the elastic moduli is as of an ideal 



deal gas, 
•i), 1), 



given by: ur — 0, and n fc — N^kT/V [see Eqs. 
( |Io| ) and (11)]. The vanishing of the of the shear modu- 
lus simply indicates that the collection of free clusters is 
a fluid. The positive contribution of the free clusters to 
the bulk modulus is crucial for the stability of the system. 
Two-dimensional Gaussian networks are stabilized in the 
presence of free clusters since K = n spc + n tc = n tc > 0. 
Three-dimensional systems are stabilized provided that 
the positive contribution of the free clusters to k over- 
comes the negative contribution of the spanning cluster. 



III. ELASTICITY OF SYSTEMS OF 
NEARLY-GAUSSIAN SPRINGS 



The elastic response of polymers and polymeric net- 
works is as of systems of Gaussian springs only in the 
first approximation. It always includes a non-linear part, 
which becomes significant when the network is suffi- 
ciently stretched, much beyond its characteristic thermal 
lengths 0Jl^] . In order to study the nature of this correc- 
tion, we consider networks of springs having the spring 
energies 



-K a f (R al3 ) 2 + -a al3 (R a,3 )\ (21) 



Our choice for the spring energy is inspired by the free 
energy of a finite long polymer chain Q , where the lead- 
ing correction to the linear relation between the force and 
the chain end-to-end vector / = KR, is a term propor- 
tional to R 2 R. The elastic energy of the system is given, 
again, as the sum of all springs energies 



<a/3> 



E 

<a/3> 
E ^ 



~1 

2 
E, 



A " ' (R ap Y + -a af3 (R a ^ 



(22) 
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We assume that E\ <C Eq , and treat the quartic term per- 
turbatively. In fact, we will make a more restrictive as- 
sumption that for each bond a af} (R al3 ) 4 < K a P (R al3 ) 2 . 
Since the quadratic term Eq does not make any contri- 
bution to the elastic constants, we will mainly focus on 
the contribution of the perturbation term, E\, to them. 

In the lowest order of a perturbation theory, the elas- 
tic constants of the network are temperature independent. 
Substituting the pair potential ( p2] ) into expression (]?]) 
for the stress tensor, and expanding this expression to 
the first order in a Q/3 , yields 





K af3 RfRf 



SE 1 



(23) 



■■0) 



where SA = A — (A)q denotes a thermal fluctuation 
of the quantity A, and ( }o denotes a thermal aver- 
age with the (unperturbed) Gaussian Boltzmann weight 
exp(— Ea/kT). er^ is the stress tensor of the correspond- 
ing Gaussian network (where a af3 = 0), given by Eq.(||), 
which can be also expressed by its value at T = 



(a/3) 



(24) 



In the above expression (^Rq^ is the ith Cartesian com- 
ponent of the bond vector R^ , at the ground state of 
the unperturbed Gaussian network. 

The next step is to substitute the pair potential j2^ ) 
into the expression for the elastic constants |l2| (see also 
Eq.(7) in Ref. [^5|). By expanding this expression to 
the first order in {a a °}, and using the fact that for the 
Gaussian network Cijki = (|20|), we find that 



ijkl 




af3 j^ap j^ap j^ap 
j k I 



Wo: 



where X is combination of terms, each of which includes 
the thermal fluctuations of some quantities. Since at 
T = there are no thermal fluctuations, that term van- 
ishes and we readily find that 



C m (T = 0) = ± ( £ a^RfRfRfR 



V ^ 



(a/3) 



PlRf 



pa/3 pa/5 p a/3 p a/3 



(a/3) 



(25) 



The second equality in the above equation is obtained by 
equating the expression inside ( )q to its value at equilib- 
rium (at zero temperature the thermal average coincides 
with this value). 



At a finite temperature we may write the elastic con- 
stants as the product of Cijki(T = 0), and a dimension- 
less function, which may depend only on terms of the 
form (kTa al3 )/{K^ s K<). Expanding the function into 
power series in these variables yields 



C m (T = 0) 



{kT a a P 
; 7 
K~< S K< 



Since Cijki (T = 0) is a linear function in the quantities 
a a @ , and since we are interested only in the first order 
correction due to the perturbation (namely, in terms lin- 
ear in a a @), we conclude that to the lowest order in a a ^ , 
Cijki are temperature independent, and therefore given 
by the above expression fl2q). 



IV. ELASTICITY OF PHANTOM PERCOLATION 
NETWORKS 

A. The percolation model 

One of the models which has been proposed to de- 
scribe the process of gelation is percolation H. In the 
percolation model, the sites or the bonds of a lattice 
are randomly occupied by, respectively, atoms or bonds, 
with an occupation probability p. In the site percolation 
model, one links every two neighboring occupied sites, 
while in the bond percolation model one assumes that all 
the sites are occupied by atoms and each pair of neighbors 
is linked if the bond between the atoms exists. Within 
the percolation model, the gel point is identified with 
the percolation threshold, the critical site/bond concen- 
tration above which a spanning cluster is formed. The 
percolation model predicts that close to the percolation 
threshold, p c , quantities like the mean cluster mass, typ- 
ical cluster linear size and gel fraction, have power-law 
dependence on (p— p c ). The relevant exponents are uni- 
versal and depend only on the dimensionality of the sys- 
tem, but not on the atomic-scale features of the system. 
The values of these exponents have been measured ex- 
perimentally for various gel systems JT6|. A fairly good 
agreement have been found between the measured expo- 
nents and their values as predicted by the percolation 
model, what proves the applicability of the percolation 
model to gelation. 

The situation concerning the elastic behavior of gels 
is not that clear. The main question is whether the 
shear modulus also follows a scaling law fi ~ (p — p c )f 
with a universal exponent /. Experimental values of this 
exponent measured for different polymeric systems are 
very scattered p7| . On the theoretical side, it has been 
demonstrated that at T = 0, the elastic behavior of per- 
colation systems depends on the nature of the interac- 
tions in the system. For non-stressed central force net- 
works the rigidity threshold occurs at a concentration of 
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bonds much larger than p c jL8| . If bond bending forces 
are present, rigidity and percolation thresholds coincide; 
however the rigidity exponent /, is considerably larger 
than the conductivity exponent, t, suggesting that the 
two problems belong to different universality classes [^9) . 
As the number of models of elasticity of random sys- 
tems increased, it became clear that de Gennes' conjec- 
ture about the identity of the exponent / to the conduc- 
tivity exponent t § can be justified only within models 
which "reduce" the thermodynamic behavior of gels to so 
called "scalar elasticity" models |^0j . Recently, the equal- 
ity f = t was measured by Plischke et al. in a numerical 
study of phantom central force percolation networks at 
T 7^ |2l]] . The authors attributed this elastic behavior 
to the entropic part of the elastic free energy. 



B. Elasticity of percolation networks 

We would like to apply our results from sections ||and 
[II to phantom percolation networks of identical springs 
having the energy E = \KR 2 (Gaussian network) or 
E 



jaR 4 (nearly-Gaussian network). We dis- 
cuss the critical elastic behavior of such networks, in the 
regime where the correlation length £ ~ (j> — p c ) u is 
much larger than the characteristic atomic length scale 
b, but much smaller then the linear size of the system 
L. The correlation length is the length scale below which 
the spanning cluster has a fractal structure and above 
which the system is homogeneous. A quantity that fol- 
lows a power law ~ (p — p c ) Y ~ when L £, 
scales as to when £ 3> L. (At p c the latter 
power law is always relevant because £ is infinite.) Since 
£ 3> b, we expect the structure of the spanning cluster 
to "forget" the details of the lattice, and have the elas- 
tic properties of an isotropic system. In the Gaussian 



case, the tensorial equality Oi 



(llST) becomes a 



scalar equality — P = S. Also, because of the vanishing 
of the elastic constants of Gaussian networks (|2^), we 
have for the shear modulus of the spanning cluster that 
/i = C44 — P = —P = £ (Q). Close to the percolation 
threshold, the conductivity scales as £ ~ (p — p c Y, and 
therefore we conclude that for Gaussian networks 



n = -p = v~(p-p c ) t , 



(26) 



in accordance with de Gennes' argument. This result is 
not changed if we also include the finite clusters, since the 
latter make no contribution to the shear modulus (just as 
they do not contribute to the conductivity of the system). 
The equality of the shear modulus and the stress, a sig- 
nature of Gaussian elasticity, was observed numerically 
in Ref. @. 

In the nearly-Gaussian case, we have from Eq.(^3|) that 
the leading term in the expression for the stress is the 
Gaussian term, and therefore we expect to have the same 
scaling behavior as in Eq. (p6|) . What distinguishes non- 
Gaussian networks from purely Gaussian ones is the non- 



vanishing elastic constants of the former. For percolation 
networks it is reasonable to assume that the elastic con- 
stants also follow a power law C ~ (p — p c ) 9 '. The elas- 
tic constants of a nearly-Gaussian networks should be 
"almost" zero, namely much smaller than the network 
stress. Therefore, the perturbative analysis in section [II 
would be self-consistent only if it yields that the expo- 
nent g > f. We can use expression (pq ) for the elastic 
constants to derive exact bounds on the value of the ex- 
ponent g. Consider a percolation network of linear size 
L in d dimensions at p c . An upper bound on the expo- 
nent g is obtained by including only a partial set of the 
bonds of the spanning cluster in the sum in expression 
(p5|). We take the set of singly connected bonds (SCBs), 
which arc such bonds that removal of each one of them 
disconnects the spanning cluster. Their number scales as 
L 1 /" §§. The force acting on a SCB is the total force 
applied on the surface of the system, which is propor- 
tional to PL^" 1 ) - £(-*/"+<*-!). The length to which a 
SCB is stretched, (Pscb)oj is proportional to the force, 
and therefore have the same scaling form 

(i2sc B )o~i ( - t/!/+ ^ 1) , (27) 
and consequently from Eq.(^5|) we get 

C ~ L- g ' v > L -d L i/v L i(-t/v+d-i)^ 

which yields the upper bound g < (4t — 1) — is(3d — 4). 
A lower bound for g is obtained by noting that for any 
bond other than the SCBs, (Pbond)o < (Pscb)o- That is 
because the SCBs are the only bonds which experience 
the total force acting on the system. We use this fact in 
expression (p5|) and write that 



C 



- —g/" 



bond)o] 

k bonds 



The term in braces in the above inequality is, however, 
proportional to the pressure [see Eq.( B4[)1 which scales 
like L~ l l v ' . This, together with result (|27j) , bring us to 
the lower bound g > it — 2v{d — 1). Using the known 
values of the exponents t and v we nn< ^ that in 

three-dimensions 2.48 < g < 2.6. In six-dimensions both 
bounds coincide to give g = 4. This last result reflects 
the fact that in six-dimensions essentially all the bonds 
of the network are SCBs. In two-dimensions we have the 
bounds 1.22 < g < 1.52. However, we must mention a 
special feature of the two-dimensional case which ques- 
tions the validity of the "nearly" Gaussian model. The 
model assumes that the contribution of the quartic term 
to the spring energy is small compared to the quadratic 
term [Eq.([j~ij)]. This happens only if the bond length 
satisfies 



i?bo„d « (K/a) 1/2 



(28) 



The longest bonds in the network are the bonds that in- 
side a cell of size £ d serve as SCBs. Close to p c , their 
length scales like 



7 



r t, V+{d -X)^ {p _ pc)t -v { <t-l) = {p _ pcf 



In two-dimensions the exponent y < 0, what implies that 
the length of the SCBs diverges, and certainly does not 
satisfy criterion (|2^). The problem is not limited to the 
SCBs only, but is relevant to a larger fraction of the bonds 
including the doubly-connected bonds, triply-connected 
bonds, and so on. It is difficult to predict, a priori, 
whether this observation should mo dify the results of the 
nearly-Gaussian model from section III. Note that we do 
not encounter such a problem for dimensionality larger 
than two, where the exponent y is positive. 



V. SUMMARY AND DISCUSSION 

We have studied the elastic properties of phantom 
Gaussian and nearly-Gaussian networks. For Gaussian 
networks, the stress and elastic constants were calculated 
exactly. We found that a characteristic feature of Gaus- 
sian networks is the vanishing of their elastic constants. 
This feature is both temperature and network-topology 
independent. We also proved the equality between the 
stress tensor of a Gaussian elastic network to the con- 
ductivity tensor of a resistor network, in which the con- 
ductance of each resistor is equal to the corresponding 
spring constant K a ^ . This result quantifies the some- 
what vague statement about an analogy between elas- 
ticity of Gaussian networks to conductivity of resistors 
networks. 

We have investigated the non-linear correction to the 
elastic behavior, by studying the properties of networks 
of springs whose energies include small quartic terms, 
additional to the leading quadratic (Gaussian) terms. 
While the stress tensor is still dominated by the contribu- 
tion of the quadratic term, the elastic constants (which 
vanish in the Gaussian network) are solely due to the 
non-Gaussian correction. We calculated the elastic con- 
stants to the first order in perturbation theory. 

Finally, we applied the results of both the Gaussian 
and the nearly-Gaussian models, to describe the elas- 
tic behavior of phantom percolation networks close to 
the percolation threshold. Obviously, the well known re- 
sult that the shear modulus follows the same scaling law, 
fi ~ (p — p c ) , like the conductivity, was recovered. We 
made a new prediction that the clastic constants also fol- 
low a scaling law C ~ (p — p c ) s ', with exponent g > t, 
and found bounds on the values of the exponent g. 

This work was supported by the Israel Science Foun- 
dation through Grant No. 177/99. 



APPENDIX A: THE CONDUCTIVITY TENSOR 
OF FINITE RESISTOR NETWORKS. 

We consider a network whose bonds are resistors of 
conductance K af3 , where the superscripts a and (3 label 
the nodes which the particular resistor connects. The 



network is finite and has an arbitrary topology, i.e., we 
make no assumption on the symmetry. We denote by B? 
the position of the node (3 and by <pP the electric poten- 
tial at the node. The network is placed inside a rectan- 
gular box of volume V = L\ x L^ x . . . x Ld, where L, is 
the length of the box along the ith Cartesian direction. 
(The derivation presented here can be easily generalized 
to systems of arbitrary shape.) The nodes of the network 
which are located on the surface of the system are called 
surface nodes, and we label them with the superscript s. 
The rest of the nodes are called the internal nodes, which 
we denote with the superscript a. The superscripts f3 and 
7 will be used to denote nodes of both types. 

The conductivity of an electrical system is a tensor, 
Sy , defined by 



(ji) = E*/ (£,-}, 



(Al) 



where the subscripts denote Cartesian coordinates and 
summation over repeated indices is implied, while (j) and 
(E) are the volume averages of the current density and 
the electric field, respectively. This definition of ap- 
plies to continuous electrical systems. It can be general- 
ized to discrete networks, if we define the current density 
by a set of Dirac (^-functions representing the currents 
in the bonds. Let us assume now that the electric po- 
tential, tp, applied on the surface of the network is such 
that on each surface point it is equal to the jth Cartesian 
coordinate of the point. Since E = — V<^, we have 



^ = vI EidV =-J& dV 



1 

v 



(pdS - 



(pdS 



where the surface integration is over the boundaries 
Xi — and Xi = Li, normal to the ith direction. How- 
ever, with our choice for the electric potential on the 



boundaries, 



j, it is easy to see that (Ei) = —Si 



where Sij is the Kronecker delta. 



The mean current density (ji), is given by 



(ji) = J Ji 



dV 



(A2) 



As we have already noted, the above definition ( [A4 ) ap- 
plies to continuous electrical systems. To make it appli- 
cable to resistor networks we need to write the current 
density as a sum of Dirac (5-functions representing the 
currents in the "linear" resistors. With this formal rep- 
resentation, the contribution to (ji) of each resistor is 
given by the line-integral 



I a>3 dx - K<*P 



R" 



bS - m 



where = K a P {ip a — ip 13 ) is the current across the re- 
sistor between nodes a and /3. Adding the contributions 
of all the resistors wc find that 



8 



(a/3) 

We may write the last result is a slightly different way 



(7i> 



^(eE* 7/507/3 



(^-^) (-J27) 



EE X7/5Q7/5 



7 



E (-«?: 



where the variable 7/3 takes the value 1 if the nodes 7 
and j3 are connected by a resistor and if at least one of 
them is an internal node; and the value 0, otherwise. The 
sums in square brackets corresponding to internal nodes 
7 = a vanish due to the Kirchoff "junction rule" for 
the vanishing of the sum of currents entering an internal 
node: 

Y^K aP Q aP (v a -yP) = 0. 


We are left with the contribution of the surface nodes 
7 = s only, i.e., 



(3i> = 



[ 8 



This last result can be also represented by summation 
over all the resistors (as), between surface and internal 
nodes 



O'i) = y 



(as) 



Finally, since the electric field is equal to (JEj) = 
we have from Eq.(Al) that 



1 

V 



E K as m (<p s - v a ) 



(as) 



We have obtained expression fll9|), which we constructed 
by mapping expression ( [l5| ) for cr^- into the electrostatic 
problem. This proves that indeed 07,- — Sy. Note that 
Sy does not depend on the positions of the internal nodes 
but only on the details of th e co nductivity. In large 
random networks the relation ( |Al|) suffices to define Sy 
without need of a detailed specification of boundary con- 
ditions. However, out exact result is valid also for small 
networks of arbitrary topology, provided that the elec- 
tric field E is generated using the very specific boundary 
conditions specified in the Appendix. 
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